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The environment of formation as a second parameter for globular 
cluster classification 
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ABSTRACT 

We perform an evolutionary multivariate analysis of a sample of 54 Galactic globular clusters 
with high-quality colour-magnitude diagrams and well-determined ages. The four parameters 
adopted for the analysis are: metallicity, age, maximum temperature on the horizontal branch, 
and absolute V magnitude. Our cladistic analysis breaks the sample into three novel groups. 
An a posteriori kinematical analysis puts groups 1 and 2 in the halo, and group 3 in the thick 
disc. The halo and disc clusters separately follow a luminosity-metallicity relation of much 
weaker slope than galaxies. This property is used to propose a new criterion for distinguish- 
ing halo and disc clusters. A comparison of the distinct properties of the two halo groups with 
those of Galactic halo field stars indicates that the clusters of group 1 originated in the inner 
halo, while those of group 2 formed in the outer halo of the Galaxy. The inner halo clus- 
ters were presumably initially the most massive one, which allowed the formation of more 
strongly helium-enriched second generation stars, thus explaining the presence of Cepheids 
and of very hot horizontal-branch stars exclusively in this group. We thus conclude that the 
"second parameter" is linked to the environment in which globular clusters form, the inner 
halo favouring the formation of the most massive clusters which subsequently become more 
strongly self-enriched than their counterparts of the galactic outer halo and disc. 
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1 INTRODUCTION 

Globular clusters are touchstones of astrophysics. galaxies. The old- 
est of them witnessed the formation and early evolution of their host 
galaxies and of their substructures, and their study has historically 
coloured the different scenarios of galaxy formation. However, as a 
collective population in a galaxy, they present unsolved problems. 
In particular their origin is not firmly established despite the large 
amount of work devoted to the analysis of correlations among their 
observable properties. It has long been realised that part of the dif- 
ficulty arises from sheer dynamical evolution undergone by these 
objects since the time of their formation. Indeed, any star cluster 
is the subject of a long list of erosive mechanisms that operate at 
different rates depending on the clu ster's location and orbi t withi n 
the Galaxy, and on its in itial mass (Djorgovski & Mevlanl d 19941) . 
iGnedin & Ostrikerl dl997I) ). 

Related to these difficulties, the search for a "second pa- 
rameter", beyond metallicity, to explain the distribution of stars 
along the horizontal branch, has met with a limited success. IZinnl 
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d 19931) suggested that the halo globular clusters break down into 
two groups according to their horizontal-branch properties: the 
two groups have different ages, kinematics and radial distributions. 
iRev et alj ( T200l|) found that an age difference can explain di fferent 
horizo ntal-branch morphologies at a giv en me t allicity . But lDotteil 
(2008), following the seminal work by iRood ( 1973) on the im- 
pact of stellar mass loss on the horizontal-branch morphology, has 
shown that considering a-element enhancement and metallicity- 
dependent mass loss along the red-giant phase produces sim ilar ef- 
fects. On the other hand lLee. Gim & Casetti-Dinescul d2007l) have 
found that globular clusters with extended horizontal branch are 
more massive than normal clusters and are dominated by ran- 
dom motions with no correlation between kinemat ics and metallic- 
ity. M ultivariate analyses have been performed bv lFusi Pecci et all 
i 19931) who find that more concentrated clusters have bluer and 
longer horizontal-branch, and by iRecio-Blanco et al. I j2006h who 
find that more massive clusters have a horizontal branch that ex- 
tends to higher temperatures. 

We can point to several reasons why previous studies have not 
been completely successful. For example, empirically separating 
the clusters into disc and halo populations solely on the basis of one 
parameter (metallicity) cannot be satisfactory, as the environment at 
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birth must i nfluence oth er cluster parameters. This was already rec- 
ognized bv lZind ( fl985T) . who combined metallicity and kinematics 
to establish two major groups (halo versus bulge/disk). The proper- 
ties of the horizontal branch have been used, but its usual character- 
isation by the parameter HBR = (B-R)/(B+V+R) is also unsatisfac- 
tory. In this respect, the parameter T e (maxim um effective temper- 
ature on the horizontal branch) introduced byi ecio-Blanco et alj 
J2006h is a welcome innnovation as shown in the present paper. 

The classical paradigm describing globular clusters as fairly 
simple systems of coeval stars of homogeneous chemical compo- 
sition has been seriously challenged recently, and this may bring a 
crucial piece to the puzzle. One fundamental characteristic of these 
systems is t heir metallicity [ Fe/H] (identified as the so-called "first- 
parameter" dvan den Bergbldl967l) )) that is generally inferred from 
their integrated photometric colours, and that varies strongly from 
cluster to cluster; in our Gal axy, glo b ular cl usters have [Fe/H] rang- 
ing between — 2.2 and (Harris! d 19961) ). Spectroscopy reveals 
that, within individual clusters, stars present very homogeneous 
contents in Fe, but also in a- and s-elements, indicatin g that proto- 
cluster s formed from gas pre-enriched in heavy metals Uames et alj 
(2004); |Prantzo7 & Charbonnel (2006)). This is in agreement with 
the predictions of quantitative models that rule out stochastic self- 
enrichment in most globular clusters as a significant contributor to 
their heavy metals, leaving pre-enrichment as the do minant con- 
tributor to metallicity [Fe/H] jBailin and Harris! {2009)). 

However globular cluster stars exhibit extremely scattered 
light-element (Li, F, C, N, O, Na, Mg, and A l) conten t s that are not 
seen among their field counterparts (see e.g. lCarrettal J2006)). This 
points to early internal chemical evolution (i.e., self-enrichment) 
in the globular cluster driven by first-generation massive and 
fast evolving stars which polluted with their hydrogen-burning 
products (among which helium in very important quantities) the 
intracluster gas out of which second g eneration stars formed 
dDecressin. Charbonnel & Mevnetl d2007h . fprantzos & Charbonnell 

( 2006) ). Recent findings of double or even multiple stellar popu- 
lations in the colour-magnitude diagrams of several globular clus- 
ters, as well as the complexity of the horizontal branch morphology 
(namely the wide colour distribution, i.e., effective temperature, of 
the stars presently burning heliu m in their co re) constitute further 
evidence for internal evolution dPiottol ( 120091) ). All these features 
can indeed be related to the presence of a second generation of He- 
enriched stars. Importantly, the star formation history depicted by 
these features seems to vary from cluster to cluster dMilone et alj 
(2008)) in a way which is still far from being understood. How- 
ever we have now firm evidence that Galactic globular clusters have 
undergone internal chemical evolution and complex star formation 
histories during their infancy that shaped their properties and in par- 
ticular their present total mass ( Decressin, Charbonnel & Mevnet 

(2007) ). This new paradigm has opened a novel route for a better 
understanding of the origin and history of globular clusters. 

It thus appeared to us that a multivariate analysis which simul- 
taneously takes into account any cosmic variance due to evolving 
physical conditions and groups objects according to environment 
of formation would be very valuable. Cladistics provides such a 
methodology. It differs from other clustering analyses in that it fo- 
cuses on evolution within a nd between grou ps rather than on sim- 
ilarities between objects (Wiley et alj (1997)). Cladistics is very 
commonly used in evolutionary biology an d has been pioneered in 
astrophysics bv lFraix-Burnet et alj d2006al lbl) : iFraix-Burnetl (2009) 
and successfully applied to the dwarf galaxies of the local group 
dFraix-Burnet eit alj d2006cl) ). 

This paper presents a multivariate analysis based on the 



method of cladistics of a large sample of Galactic globular clusters. 
After presenting the data and the method of analysis (Sect. [2}, we 
describe the three groups found by the cladistic analysis (Sect. 13. il l, 
and discuss two important results, evidence for self-enrichment 
(Sect. [4} and a possible luminosity-metallicity relation (Sect. [5). 
We then compare the properties of the three groups with those of 
Galactic halo field stars (Sect. [3J, before proposing a scenario for 
the formation of the three groups (Sect.[7}. 



2 DATA AND METHOD OF ANALYSIS 

The c hoice of parameters is a cru cial step in any multivariate anal- 
ysis. Diorgovski & Meylan (1994) have shown that the manifold of 
Galactic globular cluster properties has a dimension larger than 4, 
but that a subset of parameters linked to morphology and dynamics 
forms a three-dimensional family. Including properties of the stel- 
lar populations (e.g. a horizontal-branch paramet er or metallicity) 
will increase the manifold by 1 or 2 dimensions dFusi Pecci et alj 
dl993h ). Finally, using a large number of photometric and struc- 
tural parameters. lRecio-Blanco et alj {2006) found that 4 eigenvec- 
tors account for 79% of the total sample variance. 

Taking advantage of this indication, we selected the following 
four parameters for analysis: relative ages, metallicity ([Fe/H]), ab- 
solute V magnitude (M v ), and maximum effective temperature (T e ) 
on the horizontal branch. The age parameter is related to the sec- 
ular evolution of the stellar populations. [Fe/H] reflects the chemi- 
cal composition of the environment when and where globular clus- 
ters formed and is the "first parameter" for the horizontal branch 
morphology M v is a structural parameter that measures the present 
total baryonic mass of the globular cluster^ Finally T e is a mea- 
sure of both the pristine chemical composition of the protoclus- 
ter ([Fe/H] being the first parameter) and of the helium enrichment 
during early internal chemical evolution, since stars with higher he- 
lium content are expec ted to reach higher effe c tive temperatures on 
the ho rizontal branch jD'Antona et alj d2002h ;lRecio-Blanc o et al.1 
(2006)). The last three quantities describe truly intrinsic properties 
of globular clusters. As the age parameter evolves in all clusters, it 
cannot be used to classify them at the same level as the other three 
parameters. We thus gave it a lower weight in the cladistic analysis 
(see Appendix A). 

We performed our analysis using the large sample (54 ob- 
jects) of Recio-Blanco et al. (2006) based on homogeneous Hubble 
Space Telescope photometry. We used the T e valu e s obta ined uni- 
formly from this database by [Recio-Blanc o et all d2006t) . as well 
as the relative a ges and M v values the y adopted (i.e., taken re- 
spectively frcorijD^Angeiretal] d2005l) and the 2003 on-line re- 
vision of lHarrisI dl996l) ). We did not include more parameters in 
the cladistic analysis as we preferred to avoid the unwanted effect 
of redundancies, which give more weight to correlated parameters. 
Note also that we did not use any kinematical information in the 
cladistic analysis. However we used other parameters a posteriori 
to characterize the different groups found by the cladistic analy- 
sis: the radial velocities and structural para meters were taken from 



the 2003 on-line revision of Harris 



from Dinescu, Girard & van Altena 



1996), the orbital parameter s 



1999); [Dinescu et alj ( Piool) : 



ICasetti-Dinescu et al.1 d2007h . In this respect, the distinct orbital 



1 The absolute magnitude M v of globula r cluste rs in the Milky Way spans 
a vast range (—1.7 < M v < -10. 2, lHarrisI 419961)). and refl ects a large mass 
range (10 3 - 10 6 Mn. lMcLaughlin & van der Marell J2005ri ). 
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properties of group 3 found a posteriori (see §3) are independent 
of the methodology. 

More recent age estimation was published by 
iMarin-Franch et all {2009) after most of this project was 
completed, but for only 35 out of 54 globular clusters of our 
sample. Using these produces an inhomogenous data set, from 
di fferent sources, rely i ng on different values of Fe/H. The ages 
of IMarin-Franch et all (2009) have been determined from colour- 
magnitude diagr ams and valu e s of F e /H which are differe n t from 
those of iRecio-Blanco et alj J2006h . IMarin-Franch et all J2009h 
themselves point out the importance of using a homogeneous 
set of Fe/H to derive the ages. Nevertheless, it is instructive to 
perform analyses using both sets and compare the results, so 
as to determine how sensitive they are to the specific choice of 
parameter values. One then has the problem of combining two sets 
of ages, and it is not obvious how this should be done. It turns 
out that 9 out of 1 1 globul ar clusters that calibrat e relative ages of 
Marin-Franch et al. (2009) are in co mmon with De Angeli et al.l 
d2005l) . and using the relative ages of iDe Angeli et alj d2005l) for 
these 9 globular clusters gives a mean value of 1.0055 instead 
of 1.00, which is fine. However, comparing the ages of all the 
globular clusters in common suggests a non-linear systematic 
effect, which should perhaps be taken into account. Furthermore, 
we cannot simply convert each set of relative ages to absolute ages 
with the zero-point of each set, because the two z ero-points are 
rat her different (11.2Gyr forlDe Angeli et alj J2005h and 12.8Gyr 
for IMarin-Franch et all d2009h ). We thus simply used the relative 
ages without any attempt at homogenizing them, and the zero-point 
of one author for all ages. This additional analysis is compared to 
the main one in Sect. 13.21 

The multivariate analysis was performed using the method of 
cladistics. In short, the method works as follows. One first builds a 
matrix with values of the four parameters for each clusters. The val- 
ues must be discretized, and the number of bins (here 10) depends 
on the resolution one wants for the analysis. One then chooses a 
cluster which represents the most unevolved state, in the present 
case the metal-poorest cluster (NGC 6934), and the software clas- 
sifies all the other clusters in order of increasing diversification of 
properties (in other words, by increasing distance in the manifold of 
parameters). Clusters which are diverging from the original cluster 
in the same d irection are put on the sam e bra nch. We refer to Ap- 
pendix A and lFraix-Burnet et al.1 ( 2006a b e); iFraix-Burnetl (2009 ) 
for more details on the principles of the method. 



3 THREE GROUPS OF CLUSTERS 
3.1 The main tree 

The main result of our analysis is presented in the form of a tree 
structure, a usual form of representation in graph theory. The prop- 
erties of the sample can be read from the structure of the tree shown 
on Fig.Q] 

The tree has been rooted with group 2 which has the lowest 
metallicity on average, and as such is supposedly made of more 
primitive (or "ancestral") material. The tree divides into three main 
branches, which define three groups with quite distinct properties. 
Age increases roughly monotonically along each branch, as ex- 
pected. There are some subbranches sharing similar values of the 
adopted parameters, such as NGC 5904, 2808, 6388 and 6441, and 
which set them apart in the four-dimensional space of parameters, 
implying that they might in fact belong to a fourth small group 




Group 2 



Group 1 



Group 3 



Figure 1. The tree resulting from the cladistic analysis. The sample break 
into three distinct branches. Group 1 is composed of inner halo clusters, 
group 2 of outer halo clusters, and group 3 of disc clusters 



with properties similar to those of group 1. Specialists will im- 
mediately remark that the latter three clusters indeed share pecu- 
liar properties, in particular a very Helium-rich stellar population 
dPumo. D'Antona & Ventura! d2008h ). We emphasize that the He- 
lium abundance is not one of the parameters included in the cladis- 
tic analysis: it thus must influence one way or another the four pa- 
rameters used in the analysis. 

The properties of the three groups are presented in Table [T] 
The first 7 rows give the characteristics of the 3 groups: number 
of clusters in each group, mean distance from the Galactic center 
R gc , height above the Galactic plane Z, metallicity [Fe/H], absolute 
magnitude (M,,), mean velocity of rotati on V ro , in the Gal actic plane 
(computed with the equations given by |Frenk& White! dl98C|) us- 
ing a velocity of the Sun of 220 km/s), radial velocity dispersion a. 
The next 6 rows describe the correlations between the four param- 
eters. JTJie_bottom_£art 1 j)fd lists o r bital paramete r s take n 
fromlDinescu. Girard & van Altenal dl999l) : lDinescu et al.1 J2003h : 
ICasetti-Dinescu et al. 1 2007 ), which were not available for all clus- 
ters of each group: successively number of clusters, period of ro- 
tation P, total energy E, eccentricity of the orbit e, apocentric dis- 
tance R a , maximum distance reached above the Galactic plane Z mav , 
inclination angle with respect to the Galactic plane V, angular mo- 
mentum L, and finally the velocity components in cylindrical coor- 
dinates: vertical velocity \W\, radial velocity n and tangential ve- 
locity 0. 

Paired t-tests showed that the differences of the means of the 
groups taken two by two is not equal to (p < 0.05) for the param- 
eters log(T e ), R gc , [Fe/H], e. This is also the case between group 1 
and group 2 and between group 1 and group 3 for M r and 0, and 
between group 2 and group 3 for Z and age. There is also evidence 
(0.05 < p < 0.1) for different means of the latter two groups for 
P, E, R a , Z max and \W\, as well as between group 1 and group 2 
for age and between group 1 and group 3 for \W\. We emphasize 
that the rotational properties of Galactic globular clusters are very 
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Table 1. Properties of the three groups of globular clusters. No age estimate 
is available for NGC 6139, 6229, 6304, 6388, 6441, 6569 and 6642. The 
middle part of the Table gives the correlations: + means a correlation, - 
means an anticorrelation, x no correlation. The orbital properties presented 
in the lower part of the Table are only available for a subset of each group. 
Numbers in brackets are rms dispersions 





Group 1 


Group 2 


Group 3 




inner halo 


outer halo 


thick disc 


Number of clusters 


25 


11 


18 


R gc (kpc) 


9.4 (7.4) 


12.9 (8.0) 


4.2 (2.9) 


Z (kpc) 


4.8 (4.5) 


8.6 (7.6) 


1.9(2.0) 


Fe/H 


-1.40 (0.35) 


-1.92(0.16) 


-0.92 (0.35) 


M v (mag) 


-8.5 (0.7) 


-7.6 (0.6) 


-7.1 (0.9) 


Vrot (km/s) 


-7. 


+46. 


+ 119. 


<x (km/s) 


120(107) 


151 (107) 


69 (74) 


Age 


9.98 (0.96)" 


11.17 (0.70) 


10.18(0.48)* 


Ase- loe(T„) 


+ 


\ 




Age- [Fe/H] 






+ 


Age- My 


X 


X 


X 


log(T e ) - [Fe/H] 


+ 


X 


+ 


M r - [Fe/H] 


+ 


+ 


+ 


My - log(T e ) 


+ 


X 


X 


Number of clusters 


12 


8 


5 


P (Myr) 


353 (212) 


391 (254) 


142 (35) 


E (10 2 fon 2 /s 2 ) 


-691 (341) 


-649 (339) 


-1027 (111) 


e 


0.63 (0.18) 


0.54 (0.18) 


0.21 (0.10) 


Ra (kpc) 


17(12) 


19 (13) 


6(1.3) 


Znax (kpc) 


7.2 (6.2) 


9.6 (8.5) 


1.5 (1.0) 


<F (deg) 


32 (12) 


38 (17) 


21 (15) 


L 


886 (802) 


941 (575) 


866 (314) 


\W\ (km/s) 


100 (78) c 


77 (40) 


34 (16) rf 


n (km/s) 


+19 (141) 


-33 (124) 


-3 (22) 


(km/s) 


+15 (121) 


+89 (129) 


+170 (33) 



' average of 21 values; * average of 15 values; 
c average of 13 values; d average of 6 values. 



uncertain, since they are derived from projected radial velocities or 
numerical simulations of orbits in a model Galaxy. 

These properties show that the first two groups belong to the 
halo population of clusters, while the third group belongs to the 
thick disc population. Hereafter, the thick-disc clusters will simply 
be called disc clusters. The average velocity of rotation of group 1 
and 2 together is V ro , = 9 km/s. Group 3 is confined to the Galactic 
plane, and has a high V TOl and low a. If we separate group 3 into 
two subgroups of equal size according to their distance from the 
Galactic center, we find that V mt is 88 km/s for the inner subgroup 
(R gc < 3 kpc) and 187 km/s for the outer subgroup. There is also 
evidence that group 3 has a shorter P, lower e, *P and Z max , and no 
radial motion, as expected from clusters that partake in the overall 
rotation of the disc. 

One cluster has certainly been misclassified. NGC 6981 is in 
group 3 although it is at Z = 9. 1 kpc and has a low velocity of ro- 
tation. Since NGC 6981 is a borderline cluster in all the Figures, 
the value of one of the four parameters may be erroneous. Indeed, 
raising M v from -7.04 to -7.27 brings it into the next bin in our 
cladistic analysis (which requires that the data be discretized into 
a limited number of bins), and running the cladistic analysis again 
moves the cluster to group 2. Moving any of the other three param- 
eters by one bin and redoing the cladistic analysis does not change 
the status of the cluster. Since we found no reason for an erroneous 
My, we left it in group 3. 



Another possible discrepant cluster is NGC 6266, which is in 
group 1, but which, according to [Dinescu et al. (2003), belongs to 
a rotationally supported system, on the basis of its kinematics (but 
without precise orbital determination), despite its low metallicity. 
However, it is on the same subbranch as (and undistinguishable 
from) NGC 7089, which d efinitely belongs to the halo, according 
to our analysis and that of lDinescu. Girard & van AltemH 19991) . In 
addition, as pointed out by the referee, an isotropic distribution of 
orbits will statistically produce one or several ones in or near the 
Galactic plane. We are thus confident that NGC 6266 belongs to 
the halo. 

We now compare the statistical properties of the two halo 
groups with those of the Galactic halo field stars. The dichotomy 
of the Galactic halo stellar population has been suspected for 
some time. The mo st quantitative study in that respect, that of 
ICarollo et alj J2007t) . clearly identifies two broadly overlapping 
structural components corresponding to an inner and an outer halo. 
Stars of the inner halo are in highly eccentric orbits, in slightly pro- 
grade rotation, and have an average metallicity of [Fe/H] = -1.6. 
The outer halo stars have a uniform distribution of eccentricities, 
are in highly retrograde orbits, and of lower metallicity [Fe/H] = - 
2.2. These properties are also among those that distinguish the two 
groups of halo clusters, and indicate that the clusters of group 1 
may have originated in the inner halo: they have higher eccentric- 
ities and metallicities than group 2, while group 2 formed in the 
outer halo: they have the largest R gc , R a , Z and Z max , lowest metal- 
licities; is positive, but not significantly so. 



3.2 Additional analysis using an inhomogeneous set of ages 

We also applied our analysis to the sample of 54 globular clusters 
with composite ages as explained in Sect. [2] It leads to roughly 
the same three groups as before, with several obvious halo clus- 
ters moving into G3 (which is in principle composed of thick-disk 
globular clusters), and two globular clusters of Gl moving into G2, 
among them NGC 6218 which has a Cepheid (and should thus be 
in Gl). In other words, the number of misclassified globular clus- 
ters rises from one (NGC 6981) to only about 6. This indicates 
that our groups are fairly robust. Of course, their detailed contours 
depend on the choice of data and we believe that the slight discrep- 
ancy is due to the inhomogeneous dat a of this additional analysis: 
the ages of iMarin-Franch et alj fioogj) have been determined from 
colour-magni tude diagrams and values of Fe/H which are different 
from those oflRecio-Blanco et al.1 J2006I) . while logTe comes from 
iRecio-Blanco et a" J2006h and is derived from diagrams using the 
corresponding set of Fe/H values. 



4 EVIDENCE IN FAVOUR OF SELF-ENRICHMENT 

Self-enrichment by a first generation of stars is frequently advo- 
cated to explain the che mical anomalies in globular clusters (e.g. 
I Prantzos & Charbonnel ( 2006)) in relation with horizontal-branch 
morphology dD'Antona & Caloil (2008)). We list below distinctive 
properties of the three groups of globular clusters pointing to such 
a process. 

The well-known correlation between metallicity and extent of 
the horizontal branch is clearly present in Fig. [2] the more metallic 
the cluster, the less extended its horizontal branch at a given metal- 
licity. Also, the more luminous (that is the mo re massive) clusters 
tend to have more extended horizonta l branc h IRecio-Blanc o et ail 
(2006); iLee, Gim & Casetti-Dinescul J2007I) ). This latter point is 
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easily understood in the self-enrichment framework, in the sense 
that more massive globular clusters retain the helium-rich ejecta of 
massive polluter stars in their deeper potential well more efficiently 
than less massive globular clusters . 

The crucial new information brought by our analysis in this 
context is that the present mass of the clusters seems to depend 
on their origin, the inner halo clusters being presently more mas- 
sive than their outer halo counterparts (Fig.O. In fact, in order to 
have extreme light-element a bundance patter ns, and in particular 
extreme O-Na anticorrelation dCarrettj0 OO6)), which is linked to 



the extent of the horizontal branch I Carretta et alj l l2007h ). the inner 
halo clusters must have been even more massive in the past, before 
they lost a huge number (96%) o f first-generation low-mass stars in 
their e arly dynamical evolution jPecressin. Charbonnel & Mevnetl 
|2003)). We thus expect the mass difference between the inner and 
outer halo globular clusters to have been even larger in the past. 

The disc and halo clusters are well separated in the age- 
metallicity diagram, shown in Fig. [3] For the halo clusters, metal- 
licity decreases with age. For the disc clusters, on the contrary, it 
marginally increases with age, if at all. The spread in age is 4 Gyr 
for the halo component and only 1.5 Gyr for the disc one. The fig- 
ure confirms that the metallicity of NGC 2808 should indeed be 
about 0.5dex lower. The two other He-rich clusters of group 1 do 
not appear on this plot because no age estimate is available for 
them, but their metallicity is indeed about 0.5dex higher than the 
highest met allicity of the rest of the ha lo p opulation. This supports 
the cla im bv lCaloi & D'Antonal(f2 008) and Prant zos"& Charbonnell 
(see also lDecressin. Charbonnel & Mevnet ( 2007)) that He- 
enrichment must be associated with the build-up of abundance 
anomalies of light elements during the phase of self-enrichment. 

The presence of multiple stellar populations in the colour- 
magnitude diagram is another evidence for self-enrichment. Unfor- 
tunately, only three of the clusters of the sample, NGC 1851, 2808 
and 63 88, all belo nging to group 1, are known to have such a fea- 
ture (JPiottoj ( 2009)). We predict that two other known such clusters, 
NGC 5139, 6656 which are both metal-poor and massive, are also 
inner halo globular clusters. 

Additional clues to the self-enrichment scenario can be gath- 
ered from the RR Lyrae and Cepheid contents of the three groups. 
Globular clusters have historically been divided into two groups 
(Oosterhof I and II) according to the properties of th eir RR Lyrae 
stars. Using the compilation of Iciement et alj J200lh . we find that 
the two Oosterhof types are equally present in group 1 and 2. More 
interestingly, we find that the distribution of periods of RR Lyrae 
stars in group 2 is more sharply peaked than the corresponding dis- 
tribution for group 1, which presents a minor secondary peak at 
a higher period. The narrower period distribution in group 2 im- 
plies a small dispersion i n mass loss along the red-giant branch 
dCaloi & D'Anton^ J2008h ). while the wider distribution in group 
1 implies several generations of stars, each with a narrow distribu- 
tion of mass loss a long the red-giant branch , and with increasing 
helium abundance |D' Anto na & Caloil fc008h ). reinforcing the ne- 
cessity of self-enrichment. 

Furthermore all the population II Cepheids are found in clus- 
ters of group 1, and none in clusters of group 2. Population II 
Cepheids result from the evolution of post-horizontal branch stars 
which start from the higher temperatures of the zero-age horizontal 
branch and move toward the asy mptotic giant bran ch or leave that 
branch on rapid blueward loops JWallersteirJ J2002l) ): this explains 
their absence in halo clusters with low Te. 

We thus reach the conclusion that the inner halo favours the 
formation of very massive clusters, which retain more easily the 
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Figure 2. Metallicity-logTe diagram with symbol size indicating visual 
magnitude. The halo clusters are represented by open symbols: squares for 
group 1 and triangles for group 2. The disc clusters (group 3) are repre- 
sented by full squares 



Figure 3. Age-metallicity diagram. Same symbols as Fig. [2] No age esti- 
mate is available for NGC 6139, 6229, 6304, 6388, 6441, 6569 and 6642. 
Metallicity decreases with age in the halo clusters while it very marginally 
increases with age in the disc clusters 



products of first-generation stars and thus become more strongly 
self-enriched, giving rise to more extended horizontal branches. 



5 A LUMINOSITY-METALLICITY RELATION FOR 
GLOBULAR CLUSTERS? 

A mass-metallicity or luminosity-metallicity relation is found 
among galaxies, but is not expected and has not been found 
in the Galactic globular cluster system jDiorgovski & Mevlanl 
(1994) ), nor in numerica l simu lations of globular cluster forma- 
tion OCravtsov & Gnedinl Jiool) ). However, it has been found, in 
the form of a "blue tilt", in colour-ma gnitude diagrams of glob - 
ular cluster systems in bright ellipticals terodie & Stradeil (2006); 
lHarris et al.N2006l) : lHarriT j2009); Mi eske et aljj2006h ). Such a re- 
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lation, where metallicity increases with Mv, is present in our sample 
(Fig.[4]l, if we consider the disc and halo clusters separately. 

We have plotted a line separating the disc and halo clusters in 
Fig. [4] which can be used as a criterion for distinguishing the two 
types of clusters, in conjunction with other criteria, since the sep- 
aration is not perfect (NGC 104 is a notable exception). The disc 
clusters are on average fainter than the halo ones by about 3 mag- 
nitudes at a given metallicity. Correlation lines were adjusted to the 
two subsamples: they have comparable slopes of -2.0 and -2.8 for 
the halo and disc clusters respectively, with admittedly low corre- 
lation coefficients of 0.53 and 0.32. These slopes are much lower 
than that of - 6.75 predicted and fo und for the dwarf galaxies of the 
Local Group dDekel & Silkl ( 119860 ). They were however probably 
much steeper when the globular clusters formed, since we expect 
some of these objects to have lost a large fraction of their initial 
mass in the self-enrichment framework. 

A remarkable property of our disc clusters is that they ex- 
tend to metallicities lower than the conventional limit of [Fe/H] 
= - 0.8 dZinnl dl985D). This i s not very surprising per se, since 
iDinescu. Girard & van Altenal dl999ri have found three metal-poor 
Galactic globular clusters with thick-disc kinematics, and in M31 
there is als o evidence for metal- poor globular cl usters with disc 
kinem atics dMorrison et alj d2004l) . although see iFusi Pecci et al.l 
J200l)). 

If we divide group 3 along the conventional limit, we find that 
the metal-rich clusters have V rot = 184 km/s whereas the other ones 
have V TOt = 71 km/s, rather low, but still significantly larger than that 
of the halo clusters. We have checked that this low mean velocity is 
not due to one cluster in particular. The two subgroups do not dis- 
tinguish themselves otherwise; in particular they have the same spa- 
tial distribution (same mean R gc and Z). Two of the low-metallicity 
disc clusters, NGC 6171 and 6362 , have orbits determined by 
IDinescu. Girard & van Attend dl999k which confirm that they do 
belong to the disc population. In fact lDinescu. Girard & van Altenal 
dl999h state that their most significant result is to have shown the 
existence of metal-poor clusters with orbits consistent with the 
thick-disc motion. Our analysis confirms this finding. 

The second important result of the present paper is that the 
disc and halo globular clusters should not be separated on the 
basis of metallicity, but rather of a multivariate analysis, taking 
into account other parameters. We propose to use the magnitude 
at a given metallicity as a rough criterion, with a limit such that 
My = -2.4 x [Fe/H] - 10.45, together with other criteria, such as 
location in the Galaxy and velocity of rotation. 



6 COMPARISON WITH OTHER STUDIES AND WITH 
HALO FIELD STARS 

Before interpreting the differences between the three groups in 
terms of formation history, we com pare them to the traditional disc, 
young halo and old halo groups of lzinnl fl993h . to emphasize that 
they are rather different. There are 4, 4 and 1 young halo clusters in 
group 1, 2 and 3 respectively, and the numbers are 17, 10 and 8 for 
the old halo clusters. Two of Zinn's disc clusters (NGC 6388 and 
6441) are in our group 1. Since they are located near the Galactic 
centre and have no determined orbits, it is not possible to decide 
if they kinematically belong to the disc or the halo. Furthermore, 
as mentioned in Sect. 13.11 these two clusters might in fact belong 
to a fourth small group, maybe of Galactic bulge clusters. Turning 
to the metallicity vs HBR diagram, we find that group 1 extends 
to lower metallicities at a given HBR than the old halo clusters, 
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Figure 4. Metallicity- M v diagram. Same symbols as Fig. [2] The solid line 
is the limit between disc and halo clusters according to our criterion (M v = 
-2.4xFe/H - 10.45). The disc cluster above this line is NGC 104 

and that there are clusters of group 1 and 2 among the old halo 
clusters with the reddest HBR. Co mparing our groups to those of 
iLee. Gim & Casetti-Dinescil d2007l) . we find that all the clusters of 
their group with extended horizontal branch are in our group 1, ex- 
cept NGC 4833 (which is in our group 2). 

If we now compare (Table[2} our grouping with that of lHarrisI 

d200ll) (see his Table 1.6), we find that G3 dominates in metal-rich 
clusters class (MRC) and Gl and G2 dominate in metal-poor clus- 
ter class (MPC). In MPC alone there is no clear separation in K gc 
between Gl and G2, while G3 tends to be in the inner regions. G2 
tends to be among the more metal-poor globular clusters and tends 
to be at larger galactocentric distances. In summary, the MPC/MRC 
dichotomy corresponds roughly to our halo/disk separation, and our 
G2 populates the very metal-poor and distant MPC. Table [2] also 
shows that a classification based on arbitrary criteria does not quite 
retrieve the groups obtained with a multivariate analysis. 

We now compare the content in a-elements of the two groups 
of halo globular clusters with those of halo field stars. Sev- 
eral studies have shown that the field stars in the outer halo, 
identified through their kinematical or orbital parameters, tend 
to have l ower and/or more dispe r sed re l ative abundances in a - 
elements dStephens & Boesgaard d2002h : iGratton et at] d2003h ). 
which points to a difference in star formation rates of their birth en- 
vironment: the a-elements are indeed almost exclusively provided 
by core-collapse supernovae, which arise from ephemeral massive 
stars, while Fe is essentially produced in type la supernovae, which 
arise from stars with longer lifetimes. We find marginal evidence 
for a similar difference between groups 1 and 2 (see Fig.[5]l. [a/Fe] 
decreases with increasing [Fe/H] in group 1, while it is more dis- 
persed and shows no clear trend with [Fe/H] in group 2. This con- 
firms a similar origin for the a-elements in globular clusters and 
field stars of the halo, which thus has to occur prior to protocluster 
formation. 



7 ORIGIN OF THE THREE GROUPS OF GLOBULAR 
CLUSTERS 

It is generally assumed that the Galaxy assem bled in a hierarchica l 
fashion from collapsing haloes of dark matter dBertschingej dl998l) . 
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Figure 5. Metallicity versus light-element abundance. Same symbols as 
previous figures. The line is a least-squares fit to the globular clusters of 
group 1. Small triangles are disc (green) and halo (blue) field stars. Small 
red open squares are dwarf galaxies of the Local Group (from Venn et al. 
2004) 



Table 2. Comparison between Table 1.6 of Harris (2001) and our grouping 
by number of clusters for each class. 
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the clouds reach a critical density and are under a high external 
pressure. The mas s distribution function of th e clusters is similar to 
that of the clouds jKravtsov & Gnedinl J20Q5I) ). 

It has often been proposed that some clust ers, in partic- 
ular those identified as young halo clusters by Izinnl dl993l) . 
which tend to be in the outer halo and counter-rotating, were 
formed by accretion and disruption of satellite galaxies. But the 
chemical homogeneity of the halo, as well as substantial differ- 
ences in chemical composition between field stars in the halo 
and d warf spheroidal galaxies, arg u e against the accretion sce - 
nario ( Stephens & Boesgaard j2002h ; |Pritzl. Venn & Irwir] feOOSt) ; 
lGeisleretalJ ( l2007h ). 

The properties of our three groups of clusters can be inter- 
preted in the following way, without resorting to an external origin 
for any of the groups: 

- The clusters of the outer halo (group 2) formed during the 
initial dissipationless collapse of the proto-galaxy, from material al- 
ready polluted by earlier generations of stars, but not well homoge- 
nized and thus inhomogeneous in ff-elements. Contrary to the outer 
halo stars, they lost their initial average retrograde rotation by dy- 
namical friction and gravitational encounters. As suggested by the 
referee, this group could also have originated in some "pregalac- 
tic dwarfs" (i.e. metal-poor, gas-rich satellites that soon afterward 
began hierarchical merging). 

- The clusters of the inner halo (group 1) formed later, during 
the dissipational phase of Galactic collapse, which continued in the 
halo after the formation of the thick disc and its globular clusters. 
Since the formation of group 1 occurred later, the molecular clouds 
from which they formed had time to grow by accretion of smaller 
clumps. These clouds were already enriched at the same level in 
a-elements. Thanks to the strong potential well in the clusters (as 
evidenced by their high central velocity dispersions), the He-rich 
ejecta of first generation massive stars were not blown away and 
found their way into a more strongly helium-enriched second gen- 
eration of stars, favouring the production of hot horizontal-branch 
and Cepheid stars. 

- As indicated by their short range in age (1.5 Gyr, see Fig. [3}, 
the disc clusters (group 3) formed in a more rapid fashion than 
the two other groups, before many clusters of group 1 . This could 
presumably be due to the higher densities and external pressure 
in the thick disc. This group shows significant average prograde 
rotation, because the dissipational collapse of the disc conserved 
angular momentum. The metal-poor disc clusters seem to rotate 
more slowly and have larger eccentricities and inclinations than the 
metal-rich disc clusters. Since there is no significant age difference 
between the two subgroups, we assume that the metal-poorer clus- 
ters formed further away from the Galactic plane, and thus retained 
a larger vertical velocity component. 



Small proto-galactic clumps formed first, from initial small-scale 
density fluctuations, and collapsed in a dissipationless way (be- 
cause their gas was quickly consumed or blown away), or else 
they merged and grew in size to form larger clumps which spi- 
ralled toward the inner regions of the Galaxy by dynamical friction 
and experienced a dissi pational collap se. This scenario, or variants 
of it, first proposed by Searle & Zinn ( 1978) and verified by cos- 
mological simulations! Bekki & Chiba (2001)), has repeatedly been 
invoked to explain the properties of the Galactic globular clusters 
and of the stellar halo. Furthermore, cosmological simulations have 
shown that globular clusters can form in the densest regions of col- 
lapsing subhaloes of dark matter and giant molecular clouds, when 
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APPENDIX A: CLADISTICS APPLIED TO GLOBULAR 
CLUSTERS 

Multivariate clustering methods compare objects with a given mea- 
sure and then gather them according to a proximity criterion. Dis- 
tance analyses are based on the overall similarity derived from the 
values of the parameters describing the objects. The choice of the 
most adequate distance measure for the data under study is not 
unique and remains difficult to justify a priori. The way objects are 
subsequently grouped together (this is called the linkage) is also not 
uniquely defined. Cladistics uses a specific measure that is based 
on characters (a trait, a descriptor, an observable, or a property, 
that can be given at least two states characterizing the evolutionary 
stages of the object for that character) and compares objects in their 
evolutionary relationships. Here, the "distance" is an evolutionary 
cost. Groupings are then made on the basis of shared or inherited 
characteristics, and are most conveniently represented on an evolu- 
tionary tree. 

Character-based methods like cladistics are better suited to 
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the study of complex objects in evolution, even though the rel- 
ative evolutionary costs of the different characters is not easy to 
assess. Distance-based methods are generally faster and often pro- 
duce comparable results, but the overall similarity is not always ad- 
equate to compare evolving objects. In any case, one has to choose 
a multivariate method, and the resul ts are generally somewhat di f- 
ferent depending on this choice (e.g. [Buchanan & Collarcj I2OO8)). 
However, the main goal is to reveal a hidden structure in the data 
sample, and the relevance of the method is mainly provided by the 
interpretation and usefulness of the result. In the present paper, the 
use of cladistics is justified a priori by the evolutionary nature of 
globular clusters, and a posteriori by the strong astrophysical sig- 
nificance of the grouping found. 

In astrophysics, cladistics has already been applied to galaxies 
because they can be shown to follow a transmission with modifica- 
tion process when there are transformed through assembling, inter- 
nal evo lutio n, interaction, merge r or stripping dFraix-Burnet et al.l 
( 2006a b c); iFraix-Burnetl d2009h ). For each transformation event, 
stars, gas and dust are transmitted to the new object with some 
modification of their properties. For globular clusters, interactions 
and mergers are probably rare. It was previously thought that once 
they assembled, only the stellar ageing would affect their proper- 
ties. Nowadays, we have firm evidence that internal evolution can 
create another generation of stars, and clusters can lose mass. Ba- 
sically, the properties of a globular cluster strongly depend on the 
environment in which it formed (chemical composition and dynam- 
ics), and also on the internal evolution which includes at least the 
ageing of its stellar populations. To compare globular clusters, it is 
thus necessary to take into account the different stages of evolution 
of both the objects and their environments of formation. Since the 
clusters form in a very evolving environment (evolution of the Uni- 
verse and the dynamical environment of the parent galaxy), the ba- 
sic properties of different clusters are related to each other by some 
evolutionary pattern. In particular, the dust and gas from which the 
stars of the globular clusters form have been "polluted" (enriched in 
heavy elements) by more ancient stars, being field stars or belong- 
ing to other globular clusters. This results in a kind of transmission 
with modification process, which justifies a priori the use of cladis- 
tics. It must be clear that this is not a "descent with modification" in 
the sense that there is no replication. But evolution does neverthe- 
less create diversity. We are dealing with phylogeny (relationships 
between species), not with genealogy (relationships between indi- 
viduals). Since a multivariate classification of globular clusters is 
not yet available, we assume in the present work that each cluster 
represents a species that will have to be defined later on. 

As our work on galaxies has shown us, it is important to re- 
move parameters that are redundant. Since previous studies of the 
manifold of Galactic globular clusters have shown that 4 parame- 
ters are sufficient to describe their diversity, we selected 4 param- 
eters, 3 of which are intrinsic characteristics of the environment of 
formation. The fourth one, age, is particular in the sense that it does 
not inform on the conditions when the clusters formed, and is not 
discriminant for clustering because it evolves similarly for any clus- 
ter (parallel evolution). However, age is useful to rank the clusters 
within each group. Consequently, we applied to age a weight half 
that of the other three parameters. In addition, a stepmatrix was em- 
ployed to impose the irreversibility constraint on the age parameter 
(age can only increase). In contrast to multivariate distance meth- 
ods, undocumented values are not a problem in cladistics analyses. 
This is why the seven galaxies that have no age determination (see 
Fig[3]and Table[T]( have not been excluded in our work. 

In this paper, we use parsimony as the optimisation crite- 



rion. This works as follows. One first builds a matrix with val- 
ues of the four parameters for all clusters. The values for each 
parameter are discretized into 10 bins representing supposedly 
evolutionary states. Discretization of continuous variables is quite 
a c omplex problem, especially in the evolutionary contex t (see 
e.g. lGoloboff et alj j2006h : lThuillard & Fraix-Burned J2009l) ). The 
choice of the number of bins cannot be made in a simple objective 
way. Here, we took equal-width bins, and considered a compro- 
mise between an adequate sampling of continuous variables and 
the uncertainties on the measurements. The first constraint is given 
by the software (32 in this case). The second one would a pri- 
ori give a lower limit of something like total range / uncertainty, 
but Shannon's theorem would multiply this by 2. Hence, 10 bins 
would ac count for about 20% measu rement errors, which is rather 
large, but iRecio-Blanco et al. I d2006h do not provide precise error 
estimates, especially for logTe. Even so, border effects always im- 
ply that some objects could belong to a bin or its neighbour, a pro- 
cess that add some more artificial noise. The best way to avoid this 
effect is to make several analyses with different number of bins and 
check that the result does not depend on this number. We have done 
this for 3, 5, 8, 10, 12, 15 and 20, and the result is identical, to at 
most one misplaced cluster, for numbers higher than 8. For 3, no 
structure is found, and for 5 bins the groups are not well defined. 

Then, all possible arrangements of clusters on a tree are con- 
structed, and using the discretized matrix, the total number of state 
changes is computed for each tree. The most parsimonious tree 
is finally selected. If several such trees are found, then a consen- 
sus (strict or majority rule) tree is built. The whole procedure is 
computerised since the number of arrangements is here very large. 
The result is a diversification scenario that should be confronted 
to other knowledge and parameters. Maximum parsimony heuris- 
tic sea rches were performed using the PAUP*4.0blO (Swofforf 
(2003)) package. T he results were interpreted wit h the help of the 
Mesquite software (Maddison & Maddiso^ J2004h ). 

The tree presented in Fig. 1 is a majority rule consensus tree 
of 20000 trees, the strict consensus tree showing exactly the same 
three groups but with group 1 and 2 slightly less resolved. To fur- 
ther assess the robustness of the tree, it was not possible to make 
bootstrapping due to the irreversibility constraint on the age pa- 
rameter, and it would not have been very significant with only 4 
parameters. We performed other analyses using only 3 parameters, 
excluding the age. They all gave essentially the same three groups, 
but they were individually slightly less resolved, as expected. All 
these convergences yield strong confidence on the tree shown in 
Fig. 1. In the end, the most important point is the astrophysical in- 
terpretation we are able to give of the results. 

On Fig. 1, the tree is rooted with group 2 as an outgroup. This 
is not strictly necessary in cladistics, and we find here the same 
three groups whatever the root chosen or even on the unrooted tree. 
But we know that a low metallicity is an ancestral state for stars in 
general, this is why we have chosen group 2 a posteriori because it 
has a homogeneously low metallicity. 



